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Abstract. We discuss a class of linear control problems in a Hilbert space set- 
ting. This class encompasses such diverse systems as port-Hamiltonian systems, 
Maxwell's equations with boundary control or the acoustic equations with bound- 
ary control and boundary observation. The boundary control and observation acts 
on abstract boundary data spaces such that the only geometric constraint on the 
underlying domain stems from requiring a closed range constraint for the spatial 
operator part, a requirement which for the wave equation amounts to the validity 
of a Poincare-Wirtinger-type inequality. We also address the issue of conservativity 
of the control problems under consideration. 
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1 Introduction 



Finite-dimensional linear control problems are commonly discussed in the form of a differential- 
algebraic system. The first system equation links the state x taking values in M. n to the control 
or input u, which takes values in M m via matrices A, B of appropriate size in the way 

x(t) = Ax(t) + Bu(t), t G M>o- 



The latter equation is also known as state differential equation. This equation is completed by 
some initial condition for the state, i.e., x(0) = xq. In control theory one is mainly interested 
in the observation or output y, which is a R'-valued function given by the observation equation 

y(t) = Cx(t) + Du(t) t e M>o, 



for suitable matrices C and D. 

Thus, denoting the time-derivative by 8q and using the whole real line instead of M>o, which 
transforms the initial condition into a Dirac-<5-source term on the right-hand side, we arrive 
at the following system 




Systems of such general block structure have been generalized to the infinite-dimensional case. 
In this case A,B,C and D are linear operators in suitable Hilbert spaces. A solution theory 
for this problem is rather straightforward, if one assumes that A is a generator of a strongly 
continuous semi-group and the operators B,C and D are bounded. If one studies systems 
with boundary control, the assumption on B and C to be bounded has to be lifted. Hence, 
more sophisticated techniques need to be used to establish well-posedness of such systems, 
(T6l \T7\ [TJ [201 EH [221 E]. In this note, we shall present a unifying way of looking at 
control problems of this type which may make the solution theory more easily accessible. More 
precisely, we will provide evidence that linear control problems are nothing but evolutionary 
equations as studied in |12| . We will show that a large class of linear (boundary) control 
systems fits into this class. We exemplify these observation with linear boundary control 
problems studied by |21| \22\ [9l [101 ES]- It turns out that these equations can be discussed 
with only moderate smoothness assumptions on the boundary of the underlying domain, which 
is of prime importance in applications. Assumptions that ensure closedness of the range of 
particular differential operators are still necessary. 

A particular subclass of port-Hamiltonian systems ([SJ H]) can be discussed within this 
theory. As a by-product we give a possible generalization of these boundary control systems 
similar to port-Hamiltonian systems to the case of more than one spatial dimension. We 
will also address the issue of conservativity. Thereby, we show that the hypotheses on the 
structure of the material law in |15| can be weakened. We obtain a certain general energy- 
balance equality, imposing assumptions on the structure of the equation that are easily verified 
in applications. 

In Section [21 we give the functional analytic preliminaries needed to discuss evolutionary equa- 
tions in the sense of |12j . This includes the time-derivative realized as a normal, continuously 
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2 Functional-analytic Framework 



invertible operator and the notion of Sobolev-chains. 

Section [3] states the notion of abstract linear control systems defined as a subclass of particular 
evolutionary systems. We show well-posedness of the respective systems under easily verifiable 
conditions on the structure of the operators involved. In essence, this section recalls the well- 
posedness theorem of |12| including the notion of causality defined in [3|- 

Section [3] discusses a qualitative property of abstract linear control systems namely that of 
conservativity in the sense of [22J. In order to show conservativity of abstract linear control 
systems a particular structure of the operators involved and a regularizing property of the 
solution operator associated to the system is needed. The regularizing property is slightly 
stronger than the one in |15j. As a trade-off, the structural requirements on the operators 
involved are less restrictive. 

The subsequent section, Section [SJ provides a way to embed linear boundary control systems 
into abstract linear control systems. As a first illustrative example of boundary control systems 
we discuss in Subsection 15.11 the notion of port-Hamiltonian systems as introduced in [8j E] • 
In order to give higher-dimensional analogues for a particular subclass of port-Hamiltonian 
systems, we define abstract boundary data spaces (Subsection I5.2j> . The latter can and will 
be introduced in a purely operator-theoretic framework. Consequently, in applications these 
spaces may be defined without any regularity assumptions on the underlying domain. The 
main idea is to replace the classical trace spaces, which may not be defined in the general situ- 
ation of irregular boundaries, with an abstract analogue of "1-harmonic functions". Subsection 
15.31 provides the solution theory of a class of abstract linear control systems with boundary 
control and boundary observation. 

The last section, Section[6l is devoted to illustrate our previous findings. We give an alternative 
way to show the well-posedness of Maxwell's equation with boundary control similar to the one 
discussed in |19| (Subsection 16. 2p and the well-posedness of a wave equation with boundary 
control and observation generalizing the one discussed in [22 j (Subsection 16. ip . 



In this section we introduce the framework for evolutionary equations, which will be defined 
in the next section. The relevant statements of the results can be found in more detail in |14j . 
First, following [3], we define the time-derivative as a normal, boundedly invertible operator 
in a suitable L2-type space: 

Definition 2.1. For v £ R>o we denote by H Ut o(M.) the space of all square-integrable func- 
tion^ with respect to the exponentially weighted Lebesgue-measure exp(— 2vt) dt, equipped 
with the inner product given by 



Throughout we identify the equivalence classes induced by the equality almost everywhere with their repre- 



2 Functional-analytic Framework 




f(t)*g(t) exp(-2ut) dt (/, g G H V , (R)). 



sentatives. 
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Remark 2.2. From the definition of H u fi(M) we see that the operator exp(— um) : H Ut o(WL) — > 
L2(M), defined by (exp(— vm)f) (t) = exp(— vt)f(t), t £ R, is unitary. Furthermore, it is clear 
that the space Cq^R), the space of indefinitely differentiable functions with compact support 
on R, is dense in H Ui q(R). 

Definition 2.3. Let v > 0. We denote bjH d : if^R) C L 2 (R) ->■ L 2 (R) the usual weak 
derivative on L2(R), which is known to be skew-selfadjoint, i.e., d* = —d. We set 

9q,v '•= exp( — vtyi)~^~ {d -\- v) exp( — vrn) 

as the derivative operator on H Ut o(W). For convenience we will write do instead of do jV if the 
particular choice of v > is clear from the context. 

Remark 2.4. The operator d$^ v is normal with do tU = v. Moreover, since the operator 
exp(— i/m) _1 9exp(— urn) is skew-selfadjoint, we get that £ q{9o,v) an d H^o^ll < ^. To 
justify our choice of <9o,t> as the derivative we compute do tU (f> for (f> G (^^(R): 

(do,v4>) (<) = exp(z4) ((9 + i/) exp(-z/m)0) (t) 

= exp(z^t)(— uexp(— vm)4> + exp(— um)(j)' + z/exp(— um)4>){t) 

= m 

for all i G R. 

Next we need the concept of so-called Sobolev-chains. The proofs of the following assertions 
can be found in |144 Chapter 2]. 

Definition 2.5. Let H be a Hilbert space and C : D(C) C if — >• if be a densely defined, 
closed linear operator with £ g(C). For £ Z we set iZfc(C) as the completion of the domain 
D(C k ) with respect to the norm |C fe • \h- Then (iZfc(C)) fcgZ becomes a sequence of Hilbert 
spaces such that Hk(C) is continuously and densely embedded into Hk-i(C) for each k G Z. 
We call (iffc(C))fc e 2 the Sobolev-chain of C . We define 

tfoc(C) := p| # fe (C), 

Remark 2.6. For G N>i the operator 

C : if fe (C) -> tf fc _i(C) 
x I—?- Cx 

is unitary. For fc G Z<o consider the operator 

C : H^C) C £T fc (C) -»• fT fc _i(C) 
x h-> Cx. 

2 For the space of ^-functions defined on an open subset fi C R" with distributional gradient lying in L2(fi) n 
we use the notation 
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2 Functional-analytic Framework 

This operator turns out to be densely denned, isometric with dense range, hence it can be 
extended to a unitary operator (again denoted by C) C : H k {C) — > H k _i(C). 

Remark 2.7. 

(a) The Hilbert space H k (C) for k £ Z can be identified with the dual space H_ k (C*)* using 
the following unitary mapping 

U : H k (C) -> H- k (C*)* 

x^(y^ {C~ k x\C* k y) H ) . 

This allows an extension of the inner product (• | •) in H to a continuous sesqui-linearform 

(•!•): H k (C) x H_ k (C*) -> C 

in the sense of the dual pairing (H k (C), H_ k {C*)) . We will not distinguish between the inner 
product given on H and its extension to such pairings. 

(b) Let U be a Hilbert space and A : H\{C) — > U be a linear bounded operator. Then the 
dual operator A' : U* — > H\{C)* can be identified with the operator A° : U — > iJ_i(C*), by 
identifying the dual space U* with U and the space H\{C)* with i?_i(C*) according to the 
aforementioned unitary mapping. 

Example 2.8. Choosing H = H U $(M) for some v > and C = 8q we can construct the 
Sobolev-chain associated to 8q. We will use the notation H Vt k(M) ■= H k {do) for k £ Z. The 
Dirac-distribution 5 is an element of H Ut -i(M) and c^ 1 ^ = Xr >0 . 

Remark 2.9. For a densely defined closed linear operator A : D(A) C Hq — > Hi, where Hq 
and H\ are two Hilbert spaces, we can construct the Sobolev-chain to \A\ + i and \A*\ + i, 
respectively. Then A and A* can be established as bounded linear operators 

A:H k (\A\+i)^H k _i(\A*\+i) 

and 

A* :H k (\A*\+i)^H k _i(\A\+i) 

for all keZ. 



Remark 2.10. Let v > and H a Hilbert space. For a densely defined closed linear operator 
C : D{C) C H ^ H with £ f?(C) we consider the canonical extension 1#„ (r) <8> C of C to 
the space H v> o(M.)® H , where l_f/ I/0 (R) denotes the identity on H Vj q(M.). Analogously we extend 
do to the space H U ^(M.) <S> H by taking the tensorproduct do^ u (g) 1// with the identity 1// on 
-ff. We re- use the notation C and 8q for their respective extensions to the space H v q(M) ® ff. 
Then the operators <9o and C can be established as operators on iJ v _ 0O (R) ® H-^C) := 
UfcjGZ^fc W ® More Precisely, 

9 : fl"„,k(R) ® fli(C) ->• fl„,k-i(R) ® H 3 {C) 
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and 

C : H Vtk (R) ® Hj(C) -> fT^fcOR) ® i?i-i(C) 

are unitary operators for each G Z. As a matter of convenience, we will also write 
H Vfk (R, H) for all k G ZU{-oo, 00} for H„ tk (R)®H (or Ui & H u ^{M)®H or n/ 6Z -ff^(K)(8)iI) to 
stress the unitary equivalence of the tensor products of these Hilbert spaces with the respective 
space of (generalized) Hilbert-space-valued functions. 

3 Control Systems as Special Evolutionary Problems 

In Section we shall show that many linear control systems fit into the following particular 
class. 

Definition 3.1. Let H, V be Hilbert spaces, M ,Mi G L(H), J G L(V,H) and A : D(A) C 
H —7- H skew-selfadjoint. For v G K>o, we define the set 

{{x, /) G H v - X (K, ff V)\(d M + Mx + A)x = Jf} . 

The set £m ,Mi,A,J := Ll/>o £ M ,Mi,A,J is called evolutionary system. The system £m ,m u A,j 
is called well-posed if there exists fo G M>o such that for all v G R>i/ the relation 

{(/,x)|(x,/)G^ 0!MlAJ n^ > o(K ) iJeF)}c^ (R,F)©ij, ) o(R,F) 

defines a densely defined, continuous linear mapping from H Ut Q(M,V) to H u ^q(M, H). We call 
S^ Iq Mi A j solution operator (for v). 

Theorem 3.2 (|15t I12|). Let £m ,Mi A J be an evolutionary system. Assume that Mq = Mq 
and that there exists c G M>o such that 

vM + 9te Mi > c > 

for all sufficiently large v G M>o- Then £m ,m 1 ,a.j is well-posed and the corresponding solution 
operator Mi A j is causal, i.e., for all a G R we have 

X(-oo,a)( r no)Sl /IoMlA jX(-oo,a){ m o) = X(-oo,a) ( m o)SM ,A/i,A, j, 

where X(-oo,a)( m o) denotes the operator of multiplying with the cut-off function X(—oo,a)- 

The following proposition can be found in |15j . The basic fact, which is used in the proof is 
that 8qI commutes with S^ o m 1 A J f° r a well-posed evolutionary system £ Mq Mi a j, for all 
sufficiently large v G M>o- 

Proposition 3.3. Let £m ,Mi,AJ ^ e a well-posed evolutionary system. Then, for all suffi- 
ciently large v G K>o, we have that S^ Io Mi A j uniquely extends to a continuous linear operator 
from H„ jk (R, V) to H Vjk (M, H) for allk e'z. 

Remark 3.4. This proposition provides a way to model initial value problems, since initial 
conditions can be represented as a Dirac-source term, which turns out to be an element of the 
space H U ^\(M, H). 
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4 Conservative Systems 



We can now describe abstract linear control systems as particular evolutionary systems. 

Definition 3.5. An evolutionary system 8 M Mi A j is called abstract linear control system if 
there exist Hilbert spaces Hq, Hi,Y,U\, a densely defined, closed linear operator F : D(F) C 

/ —F* 0\ 

H -> Hi, B £ L(Ui,H) such that H = H © #i © Y, A = IF , V = H © U ± , and 

\ 0/ 

J = ( 1 B ). The Hilbert spaces Hq © i^i, U\ and y are called state, control and observation 
space, respectively. We also write Cm ,Mi,f,b to denote an abstract linear control system. 

Corollary 3.6. Let Cm ,Mi,f,b be an abstract linear control system. Assume that Mq is 
selfadjoint and that 

vM Q + 9te Mi > c> 

holds for all sufficiently large v £ R>o- ThenCM ,Ah,F,B is well-posed and causal. The solution 
operators uniquely extend to continuous linear operators from H v f.(R, H © Ui) to H U ^(M., H) 
for all k £ % and v £ R>q sufficiently large. 



-F* 0' 

Proof. Observing that | F | is a skew-selfadjoint operator, we are in the situation of 



Theorem 13.21 and Proposition 13.31 □ 



4 Conservative Systems 

In this section, we consider a qualitative property of solutions to particular linear evolutionary 
equations, namely that of conservativity. For this, a suitable regularizing property has to be 
additionally imposed. As a slightly modified version to the definition given in |15j . we define 
(locally) regularizing systems as follows: 

Definition 4.1. Let £m ,Mi,a,j be a well-posed evolutionary system. We say that Sm 0) MiA,J 
is (locally) regularizing if the following conditions are satisfied: 

(a) There exists U C D(A) dense in H such that for all T £ R and v £ R>o sufficiently large 

XR <T (m )P (0%M o + Mi + A)- 1 £<g)Mo-x K>o <8>Po) [U] C X m <t (mo) [#„,i (R, #)] , 
where Pq : H ^ H denotes the orthogonal projector onto Mq[H]. 

(b) There exists C £ R>o such that for all $ £ H we have for all T £ R and ^ £ R>o 
sufficiently large 

» <T (m ) (do, v M + Mi + ^)^ X (6 © M $) £ #„, (M, H) 

and 



X R<T (m ) ((9 M + M x + ^r 1 5 © M $) 



< C l^lo 

v,0,0 
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Remark 4.2. As we shall see in our discussion of regularizing evolutionary systems, it often 
suffices to study the following weaker norm on the left hand side of the estimate in (b): 

\f\e,u,-i,i ■= snp^H v ,i(,-e,e;H),\<i,\^i f) w ,o,o\ + \XR >e ( m )f\u,o- Then the modified inequality 
to impose is: for all T G M. and v G R>o sufficiently large and all e G M>o there exists C G M>o 
such that 

X K<T (mo) ((<9 M + Mi + ^l)" 1 5 <g> M $) 

We first will consider a conservation property for evolutionary systems. In the light of |22j 
this can be interpreted as a energy balance equality. 

Theorem 4.3. Let £m M%,A J be a regularizing well-posed evolutionary system. Let uq G H 
and consider the solution x G H u ^i(M>, H) of the equation 

(d M + Mi + A)x = 6 ® M u . 
Then the following conservation equation ZioZdjl 

/ (x\MtM lX ) = l(x\M x) (a) - hx\M x) {b) 

J[a,b] 2 2 

for almost every a, b G M>o with b > a. 

Proof. Let vo G U. Since £mq,Mi,A,J is well-posed there is a solution y G H v x($L, H) of 

(d M + M 1 +A)y = 5® M v . 

This can be re-written as 

d M (y -xr >0 ®v ) + M^y -xr >0 ®v ) + A(y -xr >0 ®v ) = -Xr >0 ®Miv -xm >0 ®Av (1) 

from which we read off that y — Xr >0 <8> & H U ^(M, H) and hence y G H^o(K,H). Let 
4> G C' 00 (]R>o) and set T := supsuppc/>. By assumption we have that XM. <T ( m o)Po(y — Xm. >0 ® 
v o) £ XlR<T( m o)[^,i(^) H)] an d hence we get from ([1]) that 

y-XM>o® w o G XR< T ( m o)[-£^,o(K,#i(A + 1))]. 

Since ^ £[/CD(i) we obtain that y G XM<T( m o)[-f^,o(^> We apply %Rt((f)y\-) Vt o,o 
to dl]) and obtain 

<Rz((f>y\doMo{y - xr >0 ® «d))i/,o,o + £Ke(&/|Afiy)„ )0j o + ^(^1^)^,0 = 0. 
Since ?/ takes values in the domain of A and since A is skew-selfadjoint, we get 

9\e(0y\d o M o {y - X R >0 ® «d))*,o,o + ^(^IM^^o = 0. (2) 
Since this holds for every eft G C' 00 (lR>o) it follows that 

<Rt(y\d M (y - X R >0 ® «fl)>o = " <He(y|M iy > a.e. on R >0 . (3) 
3 Note that Xn <T x £ .ff„,o(lR, H) for each T £ R according to the second assumption for regularizing systems. 



< cm 



E.V.-1A 
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4 Conservative Systems 



Let a, b G R >0 with a < 6. From XR< b 0™o)Fb(?/ - XR> ® ^o) G XR <6 (^o)[^,i(]R, F)] we get 
that (P y)' G L 2 ((a,b),H) with 

(-Po2/)' = <9 -Po(y - XR>o ® «o) on (a, b) 
and thus, integrating equation ([3]) over [a, 6] gives 

l(y\M y) (a) = [ (y\ <RtM iy } + Uy\M Q y) Q {b). 

1 J[a,b] 1 

Let now uq £ H and (v n ) n a sequence in U converging to no in H. For n G N let y n := 
(<9 M + Mi + A) _1 <5 ® M v„ and x := (d M + M 1 + ^) _1 (5 (g) M u . Then for every TeR 
we can estimate: 

\xr< t ( x - Un)\u,o,o = \XR< T (do M o + M i + A )~ 1 ( 5 ® M u - (5 <g> M v n )\ ufifi 
< C\u - v n \ 

where C is a some positive constant. As n — > oo we may assume y n — >• x almost every- 
where on M<b by re-using the notation for a suitable subsequence of (y n )neN an d consequently 
J^ ab Jy n \*Jit M±y n )o — > Jj a b ^(x\ 9\z M\x)q for all a,b G R. Thus, the conservation equation 
holds almost everywhere. □ 



On the Structure of Conservative Control Systems 

For the particular case of an abstract linear control systems, we shall derive now a different 
conservation property based on our observation concerning evolutionary systems. Following 
the block structure of the operator matrix A for the operators Mq and Mi we shall denote 
the corresponding entries of Mq and M\ as Mq^j and M^.,-, respectively, for i,j G {0,1,2}. 
We also have to assume a smoothing property for the solution operator. 

Theorem 4.4. Let Cmq,Mi,F-B be an abstract linear control system. Assume that Mq is selfad- 
joint and that there exists c > such that for all v > large enough, we have vMq + IHe Ml > 
c. Moreover, assume that Cm ,Mi,f,b is a locally regularizing evolutionary system and that 
Mo^o = 0, Mo 21 = 0, Mo 5 22 = 0. Assume the compatibility conditions^ 

M l 2 1 2 M i,2oV B 2 = Bo and (m^ 2 M 1i21 Y B 2 = B v 



Then for (v, w, y) G ifj, _i(R, Hq © H\ © Y) and u G H U ^(M, U±) satisfying 

—F* 0\ \ / v\ fv 
d M + Mi 4 | F ]] \ w = S ® M \ w Q | F Bu 
/ / \ W \yo 



4 Note that the condition vMq + 9\e Mi > c together with Mo, 20 = 0, Mo,2i = 0, Mo, 22 = implies that Mi, 22 
is continuously invertible. 
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for some (vq, wo,yo) E Hq © Hi © Y the control conservation equation holds: 




for a.e. a,b E M>o with a < b. 



Proof. Similarly to the proof of the conservation equation for evolutionary systems, we show 
the conservation equation stated here for initial data (vq,Wo,uo) E U, where U is chosen 
according to the definition of regularizing systems. Hence, analogously to the proof of Theorem 

—F* 0" 

]3] we get that (v, w, y) takes values in the domain of I F I and that Mq I w 







y 



locally differentiable in L/2,loc QR>o,.ff) • Let <j> E C 00 (IR > o). Then, we obtain, similarly to ([2]), 
the equation 



d M 













w 






w 











u,0,0 

+ <Rc (<f>y\B 2 u) ufi>0 



and hence 




d M 





(4) 



u) + <Rz(y\B 2 u] 
o 



almost everywhere on M>o- We aim to substitute y in the mixed term on the right hand side. 
To this end, consider the last row equation of the general system 

Using that M\ i2 2 is continuously invertible due to the positive definiteness constraint on vMq + 
yit Mi , we therefore get that 

y = -Mi2 2 M lt20 V - M^22 M l,2lW + M^B 2 U. 

Thus we have 

mt(y\B 2 u) = £He (b 2 u\ - M^ 2 M h20 v - M^ 2 M lj21 w + M^ 2 B 2 u 
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4 Conservative Systems 

= V\t /b 2 u\M{ 22 B 2 u\ - 9k (b 2 u\M^ 2 Mi >2 qv + M{ 22 M lj21 



w 



The first term on the right hand side of (jU may - using the compatibility condition - be 
computed as follows 



<Ke (v\B Q u) Q + SHc (iy|Siit) 

*Ke (Mf 2 1 2 Mi i 2ou|5 2 'u\ + <He (m{ 22 Mi,2iw\B 2 u 



Hence, 



Now, integrating equation Q over [a, 6] yields 










) 


M 








U7 





for all a, 6 positive with a < b . Using an approximation argument as in the proof of Theorem 
we get the desired assertion. □ 



Example 4.5. In |15) we studied the conservation property of the following particular system: 



/ 

do 



\ 



V 





1 


(00) 

MV 






\ 

+ 



V 

/ « 

'c 

w 

y 



o) 

0' 
1 
y/2) 



\ 





1 / 



V 



/ 

+ ( 
V 

/ 

V 

where GRAB and MV are suitable operators such that BIV* = -GRAB. It was shown that 
this system is well-posed and locally regularizing. Furthermore the compatibility conditions 
of Theorem 14.41 are satisfied with 

Mi, 22 = 1, Mi, 20 = 0, Mi )2 i = ( V2 ) , 












( zW 


\ 


N 










-y/2 j 


I 


{ or: 


) 


-1 


J 









Bo = 0, B l 



-V2 r 



Bo 



1. 
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Thus, we end up with the conservation equation 

6 



\ (Ka)| 2 + |C(a)| 2 ) -\ (\v(b)\ 2 + |C(6)| 2 ) = f \w(t)\ 2 + V2Kt( W (t)\y(t)) + \y(t)\ 2 -\u(t)\ 2 dt. 

a 

From the last row we read off the equation \/2w + y = —u and thus w = — -^(y + u). If we 
plug in this representation of w we get 



- (\v(a)\ 2 + \C(a)\ 2 ) - - (\v(b)\ 2 + |C(6)| 2 ) = J -\y{t)\ 2 - -\u(t)\ 2 dt, 

a 

which is the conservation equality in |22j . 

5 Boundary Control 

We shall now consider particular types of control equations involving so-called boundary con- 
trol. It will turn out that the notion of abstract linear control systems is rich enough to 
cover particularly interesting cases of boundary control systems such as special types of port- 
Hamiltonian systems or the control system discussed in |22| . 

As a first introductory example, we consider a particular form of so-called port-Hamiltonian 
systems (cf. e.g. OH]). 

5.1 Port-Hamiltonian Systems 

The notion of port-Hamiltonian systems with boundary control and observation as discussed 
in Section 11.2] can be described as follows: Let n £ N, a, b <E R, a < b, Po,Pi G K nxn , 
n G L 00 ((a,6),I[€ nx ' 1 ), W B , W c G K nx2n . We assume the following: 

► Pi is invertible and self-adjoint, 

► for a.e. £ G [a, b], we have ~H(C) is selfadjoint and there exist m, M G M>o such that for 
a.e. £ G [a, b] we have m < T~L(C) < M, 

( Wb \ 

► Wb and Wc have full rank and is invertible. 

V W C J 

The authors of considered the problem of finding (x, y) such that for given x^ G ^((Oj b),W' 
and u : M>o —> IK n twice continuously differentiable the following equations hold 

x{t) =P l diUx{t) + P Hx(t) 

(Hx(t))(b)\ 



y(t) =W C 
x(0) =x (0) , 





'Pi 


-Pi 






1 




'Pi 


-Pi 






1 
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where d± is the distributional derivative with respect to the spatial variable. Under particular 
assumptions on the matrices involved a well-posedness result can be obtained by using Co- 
semigroup theory, see for instance [7J Theorem 13.3.2]. Our perspective to boundary control 
systems considers a particular subclass of port-Hamiltonian (boundary control) systems. This 
subclass shows the advantage that it can be generalized to an analogue of port-Hamiltonian 
systems in more than one spatial dimension. The key assumption is that Pi is unitarily equiv- 

N* , 

n. Consequently, P\d\ is 



alent to a matrix of the form 
d!N* 

replaced by 



N 



where N £K txt with 2£ 



, with suitable domain. The unknown x decomposes into (xq,xi). 

JM u\ 

Furthermore, we assume that we only control the boundary values of x\ and that the output 
is given in terms of the boundary valued of xq. We are led to study the following problem, 
which corresponds as we will see to port-Hamiltonian systems with boundary control and 
observation as considered in [7] in a pure Hilbert space setting provided our key assumptions 
are satisfied: 

Let t G N, N G K £x£ invertible, n := 2£, M G L(L 2 ((a, b), K n )) selfadjoint and strictly 
positive definite. Let Mi G L(L 2 ((a, b), K") © K 4 ^) with the restriction of 9\e Mi to a linear 
mapping in K Ae assumed to be strictly positive definite, and Bq,B\ G lfC nxn . We define the 
operators 



Ndi: ^((0,6),] 



C L 2 ((a,b),K e ) ^ L 2 ((a,b),K e ) 
f ^ Nf, 



9iiV*: ^((0,6),^) C L 2 ((a,b),K e ) -»• £T_i(|^i| + i) 

/ 1 y (N*f)' - (N*f)(b) • 5 b + (N*f)(a) • S a . 

The expression N* f(b) is well-defined by the 1-dimensional Sobolev embedding theorem and 

N*f(b) ■ 5 b : H\(a, b),K £ ) -^K,g^ (N* f(b)\g(b)) . 

\f ^ (-Nf(b),Nf(a)), in other words C = 



with its dual, we get C° : K. £ 



We define the operator C: H\(\di\ + i) 
(-N8 b ) © N5 a . Identifying K n = K e © ] 
i), (x,y)^-N*x-5 b + N*y6 a . 

We consider the following problem: Find (xo, x\, w, y) G H U) ^i(R] L 2 ((a,b), 

,K n ) and £ ,£i € L 2 ((a, b), K e ) we have 



H-!(\di\ + 



l n ) ®K 2n ) such 



that for given u G H v $ 

( ( 

do 



V 



+ 



V 



M ,io 







-C 




M ,oi 0) 
Men 

(00) 

d x N* c° 



0) 



\ 



+ Mi 



\ 


/ x \ 








f \ 




10 


= (5 © 




+ 














B u 


/ 


v y y 








\B lU J 



°This assumptions can be guaranteed for instance for the Timoshenko beam equation, the vibrating string 
equation or the one-dimensional heat equation with boundary control, [7J. It does, however, not capture 
the one-dimensional transport equation. 
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In Section \5. 31 we shall see that this type of problem is well-posed in H u -i(M., L 2 ((a, b), ¥L n ) © 
K 2n ). For convenienc^, assume that (xq, x\, w, y) G H v> q(M., L 2 ((a, b), K n ) ©K 2n ) is a solution 





-Ndi \ ( 
-C J ^0 
(0 

[ 2n . Consequently, we get that 






o / 



/ x ° \ 

xi 
w 

V v ) 



of the above system. Then, it follows that 

takes values in L 2 ((a,b),K. n ) 

(_ 9lA r* C«)( X A €H ufi (R,L 2 ((a,b),K e )). 

Thus, with io = (wi,w 2 ) 

-N*x[ + {N*xi)(b) ■ S b - (JV*xi)(o) • <5 a - + iV*™ 2 e5 a G iJ V) o(K, L 2 ((a, b),K e )). 

The latter, however, can only happen if xi(6) = toi and x\(a) = w 2 . Hence, the first two 
equations read as 



da 



+ 



Mq oo Md m 







M ,io M ii 



+ 



-dxN* C° 
-Nd! 



Mi i00 Mi,oi M 1;02 Mi !03 
Mi,io Mi,n Mi jl2 Mi,i3 

e 



V y / 



6 



or 



/ M ,oo M 0) oi 
M ,io M 0) ii 



x 
x\ 



+ 



Mi i00 Mi, i Mi, 2 Mi )03 
Mi.io Mi,n Mi,i2 Mi,i 3 



/ x \ 

It! 

V y / 



+ 



-N*x'j 
-Nx' 



As a matter of convenience let us assume that Pi is not just unitary equivalent to 



N* 
N 



but that Pi is actually already of the form 



N* 
N 



We arrive at the following system 



M ,oo M ,oi 
M ,io Mo,ii 



x\ 



+ 



Mi i00 Mi,oi Afi )02 Mi, 03 

Mi,i M^ii Mi,i2 Ml,l 3 



Xl 
It! 

V y J 



AT* 
N 



ft 



xi 



In order to reproduce the formal structure of port-Hamiltonian systems, we are led to assume 



that 



M ,oo M ,oi 
M ,io M ,ii 



T~L and 



Mi j0 o Mi )0 i 
Mi io Mi n 



Pq. Moreover, 



Mi, 02 M 1)03 

Mi 12 Mi 13 



is as- 



^This holds true if we assume the initial data £o, £i and the control it to be smooth enough. 
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sumed to be 0. To simplify matters further, we consider the second two rows of M\ to be of 
the form 

' Mi )22 Mi )23 ' 
Afi j32 Mi )33 , 

Then the second two rows of the above system are 

Mi i22 u> + Mi j23 y - Cx = -B u 
Mi j32 u; + Mi j33 y = B x u. 

Using the above condition that ( Wl J = ( Xl (^| ] , we get that 



w 2 7 \xi(a) 



We assume that the boundary values of xq are expressed as a linear combination of the output 

Nx (b)xK >0 
-Nx (a)xR >0 

Moreover, assuming suitable invertibility properties on the operators Bq,Bi, and Mi j2 3, 

we may express the above two equations as a system of two equations of the form: 



y. Thus, there is a linear operator W € L(H Ui q(M,W 1 )) such that Wy 



u=(B - (M lj23 + W)M^B l )- 1 (M 1)22 - (Mi i23 + W)M^ 3 M 1;32 ) ) > 

y = {B B^ l M 1>m - (Mi, 23 + WO)" 1 (Mi, 22 - B B^M lj32 ) ( Xl ^X R >° ) . 

These equations are the control and the observation equations and they are of the same form 
as considered in [7] . A similar reasoning is applied in Remark 15.51 where a more general 
situation is considered. 



The discussion of boundary control within the context of port-Hamiltonian systems becomes 
accessible due to the Sobolev-embedding theorem yielding a continuous boundary trace op- 
erator and a finite-dimensional boundary trace space. In higher-dimensional situations the 
Sobolev-embedding theorem depends on the geometry of the underlying domain. A contin- 
uous boundary trace operator can only be defined for domains satisfying some regularity 
assumptions at the boundary, e.g. assuming a Lipschitz-continuous boundary. We shall ap- 
proach boundary control systems from a more general perspective without assuming undue 
regularity of the boundary. In order to have the functional analytic notions at hand to re- 
place the boundary trace space by an appropriate alias that captures the boundary data, we 
implement the necessary concepts in the next section. 
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5.2 Boundary Data Spaces 



Throughout this section, let Hq and H\ be Hilbert spaces and le10G C Hq®H\, D C H\®Hq 
be two densely defined, closed linear operators, which are assumed to be formally skew-adjoint 
linear operators, i.e. 



D C D := - ( G 
GQG :=-(£> 

Lemma 5.1. We have the orthogonal decompositions 



H 1 (\G\+i) = H 1 ( 


G 


+0 


H 1 (\D\+i)=Hi ( 


D 


+0 



[{0}] (1 - GD) 



(5) 
(6) 



Proof. Let <fi G H\ 



G 



+ i ) then for all ifi S Hi 



G 



+ i 



= fa/#)tfi(|G|+i) 

= (^>h + <|G|^||G|^ 

= m) Ho + (g^\g<p) Hi 



We read off that G<j> £ D \[G\ j = D (D) and 

DG<j> = 

The remaining case follows analogously. 



□ 



We definfH 



and obtain 



BD(G) := [{0}] (1 - DG) , 
BD(D) := [{Q}](1 -GD) 



G [BD (G)] C Si? p) , 
D [BD (25)] C BD (G) . 



7 The notation G, D is chosen as a reminder of the basic situation taking these as the closure of the classical 
operations grad and div defined on Coo— functions with compact support in an open set f2 of R n , n £ N. 
In other practical cases, these operators can change role or can be totally different operators such as curl . 

8 The notation BD ( • ) is supposed to be a reminder that in applications these spaces will serve as the spaces 
of boundary data. 
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For later purposes we also introduce the canonical projectors Kbd{G) '■ H\ ( | C| + i) —> BD (G) 
and t^bd(d) '■ H\ (\D\ + i) — > BD (D) onto the component spaces BD (G) , BD (D) according 
to the direct sum decompositions ([5]), ([6]), respectively. Note that ^* BD ^ G y 1T BD(D) are ^ e 
canonical embeddings of BD (G) in H\ (\G\ + i) and of BD (D) in Hi (\D\ + i) , respectively. 

Thus, on BD (D) we may define the operator D bj£j| 

D : BD (D) -> BD (G) 

(f>^ D(f> 

and the operator G by 

G : BD(G) -> BD{D) 
(p^Gcp. 

The operators D and G enjoy the following surprising property. 
Theorem 5.2. We have thaV^ 

G } =D = ( G X 



In particular, G and D are unitary. 

• • • • 

Proof. Obviously is DG the identity on BD (G) and GD the identity on I?D (D) . Conse- 



9 These operators are an abstract version of the Dirichlet-to-Neumann operator since the "boundary data" 
space for G is transformed into the "boundary data" space for D. Indeed, if it is a solution of the inhomo- 
geneous "Dirichlet boundary value problem" 

(1 - DG)u = 



for given data g £ BD (G) then also 

implying 
This implies 



u - g e D [G 
(l-DG) («- 5 ) = 



u = g. 



Gu = Gg 

and u is therefore also the solution of the inhomogeneous "Neumann boundary value problem" 

(1 - DG)u = 
Gu-GgeD (£>\ 

and vice versa. 
10 Note, however, that in contrast we have 

(G)* = -D 

in H fl-Dl+i) ®H (\G\+i). 
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quently, 

D = [ G 



Moreover, for (ft G BD (G) and ip £ BD (D) 

(gm) •■= (g^) =(g4>\^) +(bG(b\bA 

\ I BD(D) \ / Hi(\D\+i) \ /O \ / 



/ BD{D) \ / \ / \ / 

U\bA =: (<f>\bi/>) 

\ I ffi(|G|+i) \ / BD(G) 



leading to 

' G) =D 



in BD (D) BD (G) . □ 
Example 5.3. As an application let us calculate the dual mapping of 

7r flJ?(G) :fli(|G|+i)->5Z>(G) 

according to the Gelfand triplet H X (\G\ + i) C H (\G\ + i) C + i0 which would be 

a mapping from BD (G) (identified with BD (G) ) into H_i (|G| + i). We find 

{^BD(G)T = R* Hl (\G\+i) {^BD(G))* 
= - R Hi(|G|+i) 7r SD(G) 



|G| 2 + lW* 



BD(G) 



n BD(G) D( ^ 7r BD(G) 
^BD(G) ~ Dn BD(D)G- 



5.3 Control Systems with Boundary Control and Boundary Observation 

We apply our previous findings in this section to model problems with boundary control and 
boundary observation in more complex situations. For this purpose we consider abstract linear 
control systems Cm ,Mi,F,B where the operator F is given in the following form 

F :=(~£\ :fli(|G| + i) C ff (|G| + i) # (|£| +i) © V, (7) 

with C 6 L{H\(\G\ + i), V) for some Hilbert space V and G, D are as in Subsection 15.21 As 
a variant of |15| Lemma 5.1] we compute the adjoint of F explicitly under the additional 
constraint that G is boundedly invertible. 

"Note that the Riesz- mapping i? Hl( | G | +i) : ff_i(|G|+i) fli(|G|+i) is given by R Hl Q G \ +C) <j> — (1+|G| 2 )~V = 
(1 + G*G)"V = (1 - DG)- 1 ^. 
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5 Boundary Control 

Theorem 5.4. Let F be given as above and let G be boundedly invertible. Then 

F* : D{F*) C H (\D\ + i) © V -> H" (|G| + i) 

((,«;) ^ dc + cV 

where C° is the dual operator ofC with respect to the Gelfand-triplet H\(\G\+i) C Hq(\G\+\) C 
fl-i(|G| +i) and 

D(F*) = {(£ «/) G Hod^l + i) © ^ | Z>C + C*w G H (|G| + i)}. 



Proof. We define 



if : C flo(|I>| + i) © V -> F (|G| + i) 

(C,w) ^iJC + C^^, 

with D(A') := {(C, w) G flo(|-D| + i) © V \ D( + C°w G H Q (\G\ + i)}. From 

DC*): Hi(\D\ +i)®VC H (\D\ + i) © V -> F (|G| + i)) C A, 



we get that if is densely defined. Furthermore K is closed. Thus, it suffices to prove K* = F. 



Let v G D(K*). Then there exists 
we have 



/ 



G -Ho(|-D| + i) © V such that for all 



G D(K) 



K 



I #o(|G|+i) 

Choosing if = and £ G fli(|.D| + i) we get 

{DC\v) Ho (\G\+i) = ((\f) Ho (\D\+iy 



H (\D\+i)®V 



yielding v G ffx(|G| + i) and / = — Gv. Let now w G V be arbitrarily chosen. Like in pi 
Theorem 2.1.4] we find an element £ G Hq(\D\ + i) such that = —C^w. For this choice of 

c 



C we get (C,u>) G D(K) with if 



and thus we compute 



-Gv 
9 



H (\D\+i)®V 

(CI - Gv) Ho({f)l+i) + (w\g) v 
(D(\v) Ho (\ G \ +i) + (w\g) v 
(-C°w\v) Ho{lGl+i) + (w\g) v 
(w\ — Cv + g)v- 
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This shows g = Cv and hence K* C F. Let now v E D(F) and 



E D(K). Then 



A' 



= (DC + C te|t;)H (|G|+i) 

/ ffo(|G|+i) 



(-DCI«)flb(|G|+i) + (C°^b>H (|G|+i) 

(CI - Gv) HoQi)]+i) + (w\Cv) v , 



which shows F C K* . 



□ 



Remark 5.5. With this choice of F we can model systems with boundary observation and 
boundary control in the following way: Let Mq and M\ be of the following form 



Mn 



/M ,oo M 0j0 i 0\ 

M ,io Mo.ii 

\ 0/ 



/Mi, 00 Mi )0 i Mi, 02 Mi >03 \ 
Mi,i Mi,n Mi,i3 Mi,i 3 
M li20 Mi )21 Mi, 22 Mi >23 

\ Mi j3 o Mi, 3 i Mi )32 Mi )33 / 



for suitable bounded linear operator Mjjfc such that Mq is selfadjoint and i^Mo + M\ is 
uniformly strictly positive definite for all v > sufficiently large. Consider the abstract linear 
control system 



—F* 
d M + Mi + | F 




/ « \ 

C 

V y / 



5® M, 



Co 



+ 



\ 



o 



V ^ / 



(8) 



where F is chosen as in (J7|) and B\ E L(£/, V), S 2 E L(U,Y). We characterize the domain of 
F*. By Theorem E3] a pair (£«;) belongs to D(F*) if and only if D( + C°™ E iJod^l + i). 
Using the invertibility of D on the related Sobolev chains this is equivalent to 

C + D^C^w E H x (\b\+\). 

Hence, using the results on boundary data spaces this reads as 

K B D{D){t + b- 1 C«w) = Q. (9) 

This means that w prescribes the boundary data of £. We read off the last two lines of equation 
P and get 

M lj20 v + Mi )2 iC + M lj22 w + Mi )23 y + Cv = B x u 
Mi 30 v + Mi 3 iC + Mi i32 u; + Mi j33 y = B 2 u. 



Since the operator matrix 
assumption, we get that 



M ij22 Mi |23 

Mi i32 Mi )33 



Mi, 22 Mi, 23 
Mi 32 Mi 33 



E L(V ®Y,V © Y) is boundedly invertible by the 



Biu - (Mi j20 + C)u - Mi )2 iC 
Salt - M h30 v - M131C 
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Thus w can be expressed by v and u. If we plug this expression for w into equality ([9]) 
we obtain a boundary control equation. Likewise we may assume that the operator matrix 

i' 22 LiV © U, V © Y) is boundedly invertible and hence we get that 

-Mi,32 -D2 / 

w \ = ( -Mi )22 Si V V Mi )23 y + (Mi,2o + C)« + M lj2 iC 
it y V -Mi i32 5 2 y ^ Mi, 33 y + M li30 7; + Mi )3 iC 

This yields an expression of w in terms of y and t> and hence ([9]) becomes a boundary obser- 
vation equation. 

Example 5.6. We discuss a possible choice for the observation space, which will come in 
handy when we consider the wave equation with boundary control and observation in the 
next section. This particular choice for the control and observation space can be interpreted 
as abstract implementation of L 2 (T) of the boundary T of the underlying region. To this end, 
assume that we are given a continuous linear operator v : BD(G) — > BD(D) satisfying 



Du + u*G) 010) > (0 € BD(G) \ {0}) . 

Consider the following sesquilinear form on BD(G) : 

(■\-)u : BD(G) x BD(G) 3 (f,g) ^ ±(uf\ G g) BD(D) + \{G f\vg) BD{p y 
For / G BD(G) \ {0}, we get 

\{vf\ G f) bd(G) + \{G fWf) BD (G) = \{{bv + v* G) f\f) BD (G) > 0. 

Hence, (■\-)u 1S an inner product on BD(G). We denote by U the completion of BD(G) with 
respect to the norm induced by (-\-)u- Then U is a Hilbert space and 

j : BD(G) -> U 
/->/ 

is a dense and continuous embedding. We compute j* . Let / £ BD{G) and g S BD{G) C [/. 
Then 



(j*g\f)BD(G) = {g\jf)u 



1 • 1 • 

2^51 G f) BD(D) + 2 (G 9\vf)BD(D) 

1 • 1 • 

2^51 G f) BD(D) + 2 (" G g\f) B D(G) 

2 H G /)o + 2^ ^|/>o + G 0|/>o + g(G v* G g\G f)o 
\{{D - D)7r* BD{D) ug + (1 - DG)tt* bd ^ u* G g\f)o 

\{{D - DGD)n BD(D) ug + (1 - DG)tt bd ^u* G g\f) 
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which gives 



or 



\({l-DG)** BD{G) {bu + v* G)g\f)o 

\{{Du + u* G)g\f)o + \(G (D v + v* G)g\ G f)o 

-({Du + u* G)g\f) B D( G ), 



fg = \{bv + u* b)g 



This yields 

(D - D)ir* BD(D) G fg = \{D- t))7r BD{D) (z/+ G v* C) g. (10) 

Let us try to interpret this equation in order to underscore that this can indeed be considered 
as an equation between classical boundary traces if the boundary is sufficiently smooth. So, 
let £1 C M n be open and let grad be the weak gradient in Li2(£l) as introduced in Subsection 
16. II and let div be the weak divergence from L2(f2) n to Lzip). We denote the boundary of 17 
by r. Assume that T 7^ and that any function / € L>(grad) admits a trace /|r £ ^{T) with 
continuous trace operator. Moreover, assume that there exists a well-defined unit outward 
normal n : T — > M n being such that there exists an extension to O in a way that this extension 
(denoted by the same name) satisfies n £ Lo^fi)™ with distributional divergence lying in 
Loo(r2). Then the operator v: grad | + i) — > H\{ \ div | + i), / 1— > nf is well-defined and 
continuous. For the choices D = div, G = grad and v = 11 ' BD{Am)^' K BD{graA) m GS) we can 

interpret fjlUp as the equality of the Neumann trace of G j*g and the trace of g. Indeed, for 
f,g£ BD(gr&d) we compute formally with the help of the divergence theorem 

J grad fg ■ n f = ^ J vg ■ n f + grad g\f) + ^(gradz/* grad g\ grad/) 

r r 

If 1 

= 2 / u 9 ■ n f + -(^* gTadg\f) Hl{lgI3udl+i) 

r 

If 1 

= 2 / V 9 ' n f + div W"))ifi(| grad |+i) 

r 

= \ J U9 ' n f + \ j diy (9vf) 



r 



\jgf + \jg{vf)-n 

: 

/ 



r 

gf- 

r 
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6 Some Further Applications 

6.1 Boundary Control and Observation for Acoustic Waves 

We introduce the operator 

grad : D(grad) C L 2 {9) -»• L 2 (0) n 

as the usual weak gradient in L2^l) for a suitable domain $7 C W 1 . We require that the 
geometric properties of f2 are such that grad is injective and that the range grad[L2(0)] is 
closed^! in L2{£l) n ■ We choose to use this assumption to avoid technicalities. If grad is not 
injective, one has to proceed similarly to the way presented in the next section. However, the 
assumption on grad[L2(fi)] C L,2(£l) n to be closed is essential. See also the discussion in )18[ 
Remark 3.1(a)]. We denote by 7r gra( j : L2(£l) n —> grad[L2(^)] the canonical projector induced 
by the orthogonal decomposition of Li2{£l) n with respect to the closed subspace grad[L2(f2)] 
and consider the operator 7T gra dgrad : -D(grad) C L2(£l) — > grad[L2(fi)]. The negative adjoint 
of this operator is given by div7r* rad : D(div) n grad[L2(ri)] C grad[L2(^)] — > ^2(^)5 where 
div is defined as the closure of the divergence defined on the space of test functions C 00 (Q.) n . 
In |22l Section 7] a control system for the wave equation has been discussed, which has its 
first order pde analogue in the system: 



/ / 

do 



\ 



V 







-TTgrad grad 

c 





0) 

1 

0) 

div| 



\ 

+ 



grad[L 2 (f7)] 




0) 

/ 



V 



0) 
0' 
1 

M) 

\\ 







+ 



/ 



/ V \ 



s (0) 






+ 



V 




-V2 

-1 



(11) 



Using the Hilbert space U from Example 15.61 we define the operator C b\l^l 



C-.HtQ grad \+i)->U 

u -bjir BD ( grad )U, 

where b G L(U). Then we are in the situation of Theorem 15.41 and hence Corollary 13.61 is 
applicable. The underlying space of equation (jlip is given by H = L2(£l)©grad[L2(f2)]©?7 @U. 

12 This holds if a Poincare-Wirtinger-type inequality holds, which is for example the case, if fl is connected, 

bounded in one direction, satisfies the segment property and possesses infinite Lebesgue-measure. 
13 Note that |7r gra d grad | = | grad |. 
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We compute C° with respect to the Gelfand-triplet i?i(|grad| + i) C J^q(| grad | + i) C 
-ff_i(|grad| + i). For u G i?i (| grad | + i),v G BD(grad) C U, using Example 15.31 we get 
that 

(-C%\u) = -(v\Cu)u 

= {v\bjTi:BD(grad)U)u 

= (j* b* V \tT BD ( grad ) U) BD ( grad ) 

= ( 7r BD(grad)i* 6 ^l^)0 

= ((^BD(grad) ~ div7r BD(div) g™d)j*b*v\u) 

= ((^BDCgrad) divgrad -divvr|j D(div) grad) j*b*v\u) 

o • 

= ((div-div)vr|j D{div) grad j*b*v\u} 



and we read off that C*v = f — (div — div)^^^ grad j* b* j v G H_i(|grad| + i) for all 
v G -BD(grad) C U. Hence, using @, we write the boundary equation as 

^BD(div) U ~ d[v \grld[L 2 (Q)] (( div - div ) 7r BD(div) g rad f b *J ™) = °- 

d[v \grld[L 2 (n)] (( div - div ) 7r BD(div) S rad f b *) w e #i(| div I + i), 



Since 



we get that 

K B D(div)( = ^Bi?(div)dfypid[£ a (n)] ^(div-div)vr^ D(div) grad fb* 
= — grad j*b*w. 

To invoke the boundary control and observation equation we compute 



w 



1 0\ 1 (-y/2u-Cv\ ( 1 0\ f-V2u-Cv 



y \V2lJ I -« / \-V21 



-ii 



and 

-l 

w 



1 v 7 ^ [Cv\_f 1 -y/2\fCv 

V2 i J v y J~ \ -V2 1 J \ y 



Thus, we get w = —\^2u — Cv and w = Cv — y/2y. This yields 

n B D(div)( = V 2 grad j*b*u- grad j*b*bjir BD{gr&d) v 

and 

7i"BB(div)C =grad j*b*bjiT BD{grad) v + V2 grad j*6 y. 
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6 Some Further Applications 

Remark 6.1. Let us assume that there exists a outward unit normal n on T := Q \ Cl such 
that there exists a bounded, measurable extension to f2 with bounded, measurable distribu- 
tional divergence. Using the interpretation from Example l5.6l the assumption b*u, b*Cv, b*y € 
BDferadj^and imposing suitable additional requirements on the underlying domain, we can 
interpret the latter equations as 

n ■ C = -b*bv + V2b*u 
n-( = b*bv + V2b*y 

on r as boundary control and boundary observation equation, respectively. These correspond 
to the boundary equations originally considered in |22[ Section 7]. 

6.2 Boundary Control for Electromagnetic Waves 

As a second example we consider a boundary control problem for Maxwell's system. We shall 
first introduce the operators involved. Throughout let C ]R 3 be an open domain. 

Definition 6.2. We define the operator curl as the closure of the operator 

Coo(n) 3 ci 2 (fi) 3 ^£2(o) 3 

/ o -da d 2 \ /<M 

\-d 2 d 1 / \<fe/ 

where di denotes the partial derivative with respect to the i-th coordinate. The operator curl 
turns out to be symmetric and we set curl := ^curl^ and obtain the relation 

curl C curl . 

In |19] the exact controllability of the following problem was considered 

d eE + curl H = 5® , 
d fiH-cmlE = 5®H w , 

where the control u € i?D(curl) prescribes the boundary behaviour of the tangential compo- 
nent of H, i.e., T^BDicnr\)H = u. This problem can be dealt with in the following way: We 
introduce the function H := H — ^BD^url) 11, an< ^ formulate Maxwell's equations for the pair 
(E, H) as follows 

d eE + curlF = 5 <g> £ ,(0) - curl7r^ D ( curl) n, 
doiiH - curl£ = 5® H {0) - dofj,n% D(f:ulV) u 

14 In [22] these assumptions are formulated with the help of a certain quotient space Zq. 
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or in matrix-form 



6.2 Boundary Control for Electromagnetic Waves 



>(s:) + (-it))(i)-(5!)-(ss: 



u. 



By our general solution theory (Theorem 13 . 21 and Proposition [373]) this system is well-posed and 
we obtain a unique solution (E,H) £ H v —i(K; L 2 (JY) S © L2(^) 3 )- Since the time derivative 
of u occurs as a source term, we obtain a regularity loss of the solution (E,H), although 
the system is locally regularizing. In order to detour this regularity loss, we may follow the 
strategy of Subsection 15.31 and point out, which type of boundary control equations can be 
treated in this way. 

In the framework of Subsection 15. 3[ we want curl to play the rol^ll of D and — curl that of 
G. In view of Theorem 15.41 we have to guarantee that curl is boundedly invertible. For this 
purpose we consider the restriction of the operator curl given by 

curl : D(curl) n ([{0}] curl)" 1 C ([{0}] curl) ± -> curl[L 2 (n) 3 ]. 



We require that f2 has suitable geometric properties such that curl[L2(fJ) 3 ] is closed in order 
to obtain a boundedly invertible operator 1^*1 An easy computation shows that ^curl^ = 
curl| cur i[^ 2 (Q)3]. We decompose the Hilbert space L 2 (S1) 3 into the following orthogonal sub- 



spaces 

\3 



L 2 (ny = [{0}] curl ©([{0}] curl) 
L 2 (n) 3 = [{0}]curlffi f[{0}]curl 



and denote by tt : L 2 (Q) 3 -»• [{0}] curl, tti : L 2 (^) 3 -> ([{0}] curl)" 1 , tt : L 2 (nf -»• [{0}]curl 



and 7Ti : L 2 {^1) ~~ ^ (J{0}] curi J the respective orthogonal projections. Since T^BD{cva\)H = 
i^BD{cm\)^iH for each H £ D(curl) we may write the boundary control problem in the follow- 
ing way 



( irieirl 







tti/ittI 

(00) 
tto/j-ttI 



7ri£7r 









( \ 





V 



This implies G = —curl and D = curl . 

For example, domains ft C R 3 with conical points, wedges and cups with a cross section satisfying the 
segment property. In a large class of such domains is characterized for which the compactness of 
the embedding D(curl) (~1 D(div) <—¥ L2(0) 3 holds. This compact embedding result implies the desired 
properties for curl. 
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6 Some Further Applications 



+ 



+ 






M lj3 i 







-curl \ 
-C ) 







(00) 


Af 1>32 Mi, 33 
(00) 

I°v 

curl)* C° 











J 
V *o# (0) / 



+ 






M li34 
















Mi, 35 





/ 



\\ / TT\E \ 

w 

TT E 







o / 



where 



C : #i |curl| + i 



U 



is a bounded linear operator, U an arbitrary Hilbert space and B E L(U). The linear operators 
Mi 3i for i E {1, . . . , 5} are bounded in the respective Hilbert spaces and 9\e Mi 33 is assumed 
to be strictly positive definite. Since curl is boundedly invertible, Theorem 15.41 applies and 
Corollary 13.61 yields the well-posedness of the control problem. The domain of ^ fcurl^ C 
reads as 



TT 1 H\ 



w 



E D 



curl C c 



<^ (curl) TdH + C°w E flo(|curl| + i) 
O (curlj UriH + (^(curl 



«4> ttiH + ( ( curlj 



& ^BD(curl) ( Wl-H" + ( (^CUrlJ 



E #1 

/. N * N -1 



C°w ) E # (|curl|+i) 
+ i N 



( curl J 







(12) 



for each w E U, H E Z^fi) 3 . By the 3 rd equation of the above boundary control problem, we 
get that 



w 



-Mf^g ((Mi )31 - C)mE + Mx^iH + M h3A TT E + M h35 7r H - Bu) 



and hence (|12|) yields 
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^D(curi) UiH - ( (curl) *) C*M^ ((Mi )3 i - C)mE + M^foH 
+ M h3A ir E + M h35 Tr H -Bu\j=0 

Although this equation covers a number of possible control equations, it appears that in this 
setting the term (M±^i — C)ttiE cannot be made to vanish, since we have to assume that 
Ml 31 is bounded on i?o(|curl| +i) whereas in general C is not. This shows that in this setting 
only boundary control equations containing terms in tt\H and ttiE can be treated without 
more intricate adjustments. 
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